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ABSTRACT 

We present measurements of the ionising ultraviolet background (UVB) at z ~ 5—6 using 
the quasar proximity effect. The fifteen quasars in our sample cover the range 4.6 < z q < 6.4, 
enabling the first proximity effect measurements of the UVB at z > 5. The metagalactic 
hydrogen ionisation rate, Tbkg, was determined by modelling the combined ionisation field 
from the quasar and the UVB in the proximity zone on a pixel-by-pixel basis. The optical 
depths in the spectra were corrected for the expected effect of the quasar until the mean flux 
in the proximity region equalled that in the average Lya forest, and from this we make a 
measurement of Tbkg- A number of systematic effects were tested using synthetic spectra. 
Noise in the flux was found to be the largest source of bias at z ~ 5, while uncertainties 
in the mean transmitted Lya flux are responsible for the largest bias at z ~ 6. The impacts 
of large-scale overdensities and Lyman limit systems on Tbkg were also investigated, but 
found to be small at z > 5. We find a decline in Tbkg with redshift, from log(Tbkg) = 
-12.15 ± 0.16 at z ~ 5 to log(r bk g) = -12.84 ± 0.18 at z ~ 6 (la errors). Compared to 
UVB measurements at lower redshifts, our measurements suggest a drop of a factor of five 
in the H i photoionisation rate between z ~ 4 and z ~ 6. The decline of Tbkg appears to be 
gradual, and we find no evidence for a sudden change in the UVB at any redshift that would 
indicate a rapid change in the attenuation length of ionising photons. Combined with recent 
measurements of the evolution of the mean free path of ionising photons, our results imply 
decline in the emissivity of ionising photons by roughly a factor of two from z ~ 5 to 6, albeit 
with significant uncertainty due to the measurement errors in both Tbkg and the mean free 
path. 
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cosmology: early Universe 



1 INTRODUCTION 

The metagalactic ultraviolet background (UVB) is the net radia- 
tion field responsible for keeping the Universe ionised from the 
end of reionization to the present day. The relative contributions 
from galaxies and quasars, as well as filtering by the intergalac- 
tic medium (IGM) itself, determine the intensity and spectrum 



* A part of the observations were made at the W.M. Keck Observatory 
which is operated as a scientific partnership between the California Insti- 
tute of Technology and the University of California; it was made possible 
by the generous support of the W.M. Keck Foundation. This paper also in- 
cludes data gathered with the 6.5 meter Magellan Telescopes located at Las 
Campanas Observatory, Chile, 
f E-mail: acalver@ast.cam.ac.uk 



of the UVB (e.g. iBechtold et al.1 1 19871: lHaardt & Madaul 1 19961 : 
iFardal et alll998tlHaardt & Madaull200ll) . Thus by measuring the 
UVB one can hope to place constraints on the evolution of the 
source population with redshift. 

Of particular interest is the evolution of the UVB at z ~ 6. 
The appearance of Gunn-Peterson (GP) troughs in the spectra of the 
highest-redshift known quasars has been interpreted as evidence for 
a sharp dow nturn in the UVB at z > 6 signalling the end of reion- 
ization (e.g jFan et al.ll2006bl) . However, the diminishing transmit- 
ted flux is also co nsistent with a more slowly-evolving UVB an d 
IGM density field dBolton & Haehneltll2007bl : lBecker et aT1l2007l) . 

The bulk of the ionising photons at z > 6 that make 
up t he UVB are believed to come from low-luminosity galax- 
ies l lRichard et ail |2006t IStark et alj 120071 ; iRichard et all 120081 : 
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ISrbinovskv & Wvithel201ol ; IOesch et alj2010h . Direct searches for 

these sources at 2 > 6 — 10 have taken advantage of recent 

very deep optical and near-infrared imagine from both the ground 

H ■ 11 1 J -° — a — 1| -\ J a J 

and space JBunker et a l, 2004: Bouwens et al. 2006; Yoshi da et al.l 

120061 : iBouwens et al.ll2008l ; lOuchi et alj|2009h - The majority of 
these faint gala xies, however, still re main below current detection 
thresholds ("e.g. lBouwens et al J2010l> . 

The quasar proximity effect has been a classic tool for di- 
rectly measuring the intensity of the UVB at high redshifts. Since 
quasars are highly luminous, their output of ionising photons will 
dominate over that of the UVB out to large (up to several proper 
Mpc) distances. This produces a region of enhanced transmissio n 
near the redshift of the quasar, first noted bv lCarswell et al.l J1982D . 
known as the 'proximity region.' The size of this region depends 
both on the quasar luminosity and the intensity of the UVB. For 
a known quasar luminosity, therefore, the UVB can be estimated 
by measuring the extent of the proximity zone. Classically, the 
proximity effect has been measured by comparing column densi- 
ties of the H 1 Lyq absorption lines in the forest with those close to 
the quasar (Murdoch et alll986l:lTvtledl987l;ICarswell et al.l 19871; 



IScottetalj feoOO; Bait lik et al.lll988l . hereafter BDO), although a 
variety of flux statistics hav e also been used llLiske & WilligeJ 
120011 : iDall'Aglio et"ai] 120081 : iDall'Aglio & Gnedinl |2010|) . Alto- 
gether, proximity effect stu dies have delivered me asurements of 
the UVB from z ~ 0.5 l lKulkarni & Fail Il993h to z ~ 4.5 
( IDall'Aglio etal]|2009l) . 

More recently, an alternative method of estimating the 
UVB has been developed which uses the mean flux in the 
Lyq forest in combination with numerical simulations. The UVB 
in the simulation is adjusted until the mean flux in artifi- 
cial Lyq forest spectra is equal to that in the real data (see 
e.g.. iRauch et all Il997l: ISongaila et alj 1 19991: iTvtler et al.l 120041 : 



iBolton et alj 120051 : iKirkman et alj 120051 ; Ijena et alj 120051) . At 

z > 4.5 , the UVB has so far only been det e rmined using this 
method llMcDonald & Miralda-Escudl 1200 ll; iMeiksin & White] 
l2004l : lBolton & Haehneltll2007bl : IWvithe & BoltorJl2010l) . At 2 ~ 
5 — 6, however, converting the mean flux into an ionisation rate de- 
pends sensitively o n modelling the gas density distribution at very 
low densities (e.g. iMiralda-Escude et alj 2000h. Not only is this 
a challenge numerically jBolton & Becked 120091) . but the optical 
depth distribution will depend on the properties of the simulation, 
including the gas temperature. 

At z < 4.5, estimates of the UVB from the proximity ef- 
fect and flux decrement methods can be directly compared. These 
are generally discrepant at 2 < z < 4, with the proximity effect 
estimates systematically higher for most studies. Recent proxim- 
ity effect papers suggest that there is a competition between the 
effect of the enhanced intensity of the UVB and the overdensity 
of matter close to the quasar. Classically, the proximity effect as- 
sumes that the density distribution close to the quasar is the same 
as that in the general IGM. Not accounting for overdensities may 
therefore cause the UVB to be overestimated by up to a factor of 3 
dLoeb & Eisensteinll995h . B y using an independen t measure of the 
UVB from flux decrements, iRollinde et alj J2005T) found tentative 
evidence that quasars at z = 2 — 3 may reside in haloes as massive 
as ~ 10 14 Mfl, similar to the most massive halo in the Millennium 
simulation. iGuimaraes et ail J2007t) have claimed to detect a sim- 
ilar effect at 2 ~ 4. Other recent papers, however, have claimed 
to have overcome the environ mental bias of an enhan ced average 
density in the proximity zone. iDall' Aglio et all d2008h report that 
with their measurement method only 10 per cent of their sample at 
2 < z < 4.5 showed significant excess absorption attributable to 



an overdense environment, and that for the majority of the quasars 
in their sample their proximity effect measurement of the UVB did 
not appear to be affected by an overdense quasar environment (at 
least on scales > 3 Mpc). 

In this paper, rather than identify individual lines, which be- 
comes increasingly difficult at z > 4, we further develop a variant 
on the 'flux-transmission' method l lLiske & Willigej200lh to mea- 
sure the UVB intensity at z > 4.5. At lower redshifts this has been 
used to compare the mean flux averaged over extended sections of 
the spectrum near the quasar redshift to the mean flux in the forest. 
Rather than compute the mean flux in sections, we consider indi- 
vidual pixel optical depths across the proximity region. The optical 
depths are modified to remove the presumed effect of the quasar, 
until the proximity region has the same mean flux as the forest at 
that redshift. The characteristic scale length of the quasar model 
is then combined with the quasar luminosity to estimate the UVB. 
This simple approach avoids a direct dependence on simulations as 
it does not require the optical depth distribution of the forest to be 
known a priori. We use, however, simulations extensively to esti- 
mate the bias and uncertainties of our method. 

The remainder of this paper is organised as follows: In sec- 
tion 2 we describe the observational data as well as the hydro- 
dynamical simulations used. In section 3 we detail the proximity 
effect analysis, and the sources of systematic bias such as the ef- 
fect of the quasar environment. The results are presented and their 
implications discussed in Section 4. Finally, we present our conclu- 
sions in Section 5. Throughout this paper we assume a flat Universe 
and cosmologica l parameters taken fr om the mean of the WMAP 
5 year data set jKomatsu et alj 12009). with Hubble constant Ho 
= 72 kms -1 Mpc -1 and density parameters (fi m , Oa) = (0.26, 
0.74). 



2 DATA AND MODELS 
2.1 Observed spectra 

The quasar spectra used in this paper were taken with either the 
Keck or Magellan telescopes. The majority of the quasars at z > 5 
were ob served with the High Resolution Echelle Spectrometer 
(HIRES; IVogt et al.|[l994h on the 10 m Keck I telescope, and re- 
duced using a custom set of I DL routines and opt i mal s ky sub- 
traction techniques as detailed in iBecker et alj J2006l.l2007h . All of 
these observations were made with the 0.86" slit with R — 40 000, 
and so the velocity resolution is 6.7 km s -1 . 

The majority of the quasars at z < 5 were observed using the 
6.5 m Magellan-II Clay and the Magellan Ina mori Kyocera Echelle 
(MIKE) spectrograph l lBernstein et al.1120031) and reduced with a 
similar custom pipeline. The velocity resolution is roughly half that 
of the HIRES spectra at 13.6 km s _1 . A list of the targets is shown 
in Table[0 

Quasar redshifts we re taken either fr om the CO and 
Mg 11 redshifts presented in ICarilli et all d20ld) . or from the spec- 
tra themselves by identifying the redshift at which the Lyq forest 
appears to start. Errors on the redshifts measured from the apparent 
start of the Lyq forest were estimated by comparing to more precise 
redshifts from Mg 11 and CO where available, or those in the SDSS. 
All the objects have photometry in the SDSS, and those at z < 5.5 
also have flux-calibrated spectra in the SDSS archive. Continuum 
mag nitudes for those at z > 5.5 were tak en from the discovery pa- 
pers dFan et al J200ll, 1 200311 20041. 1 2006bh . whilst fluxes for those at 
z < 5.5 were measured from the SDSS spectra. In both cases the 
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Table 1. The list of quasars included in this paper. Columns give the quasar name and redshift, details of the observations, and the average signal-to-noise per 
pixel in the 40h~ 1 comoving Mpc closest to the quasar redshift, after masking out skyline residuals (see Section l3~3V 



Name 



Inst. 



Dates 



t eX p (bis) Ref. S/N 



SDS's T1 

O.L/iJ>J J 1 ItOTJZJ 1 


6 42 


HIRES 


Tan 9DOS - Fph 9DOS 


14 2 a 


1 


16 


snss no30+05?4 


6.31 


HIRES 


Feb 2005 


10.0 


[ 


12 


SDSS J1623+3112 


6.25 


HIRES 


Jun 2005 


12.5 


1 


11 


SDSS J0818+1722 


6.02 


HIRES 


Feb 2006 


8.3 


2 


12 


SDSS J1306+O356 


6.02 


MIKE 


Feb 2007 


6.7 


4 


14 


SDSS J0002+2550 


5.82 


HIRES 


Jan 2005 - Jul 2008 


14.2 


1.4 


21 


SDSS JO836+O054 


5.81 


HIRES 


Jan 2005 


12.5 a 


1 


19 


SDSS J0231-0728 


5.41 


HIRES 


Jan 2005 - Feb 2005 


10.0 


1 


14 


SDSS J1659+2709 


5.33 


HIRES 


Sep 2007 - Jul 2008 


11.7 


3 


32 


SDSS J0915+4924 


5.20 


HIRES 


Feb 2005 


10.0 


1 


23 


SDSS J1204-0021 


5.09 


HIRES 


Jan 2005 - Feb 2005 


6.7 


1 


17 


SDSS J001 1+1440 


4.97 


HIRES 


Sep 2007 


6.7 


3 


47 


SDSS J2225-0014 


4.89 


MIKE 


Oct 2007 


5.0 


4 


23 


SDSS J1616+0501 


4.88 


MIKE 


Mar 2008 


3.3 


4 


21 


SDSS J2147-0838 


4.59 


MIKE 


Oct 2007 


8.3 


3 


51 



a The present reductions include only data taken with the upgraded detector. 
References: 

1 - Becker et al. (2006) ; 2 - Becker et al. (2007) ; 3 - Becker et al. (2010) ; 4 - This paper 



Table 2. Mass resolution and box size (comoving) of the hydrodynami- 
cal simulations used in this work. Model C was primarily used to simulate 
quasar proximity zones. 



Model L 


Total particle 


-^gas 


[Mpc/h] 


number 


[M /h] 



A 


20 


2 x 100 3 


1.03 x 10 8 


B 


20 


2 x 200 3 


1.29 x 10 7 


C 


20 


2 x 400 3 


1.61 x 10 6 


D 


40 


2 x 200 3 


1.03 x 10 s 


E 


40 


2 x 400 3 


1.29 x 10 7 


F 


80 


2 x 400 3 


1.03 x 10 s 



continuum flux was measured at a rest wavelength of 1280 A. The 
error in the flux measured from the SDSS spectra is conservatively 
taken to be 0.5 x 10~ 17 erg cm~ 2 s _1 A -1 . 

All th e spectra we r e norm alised following the method de- 
scribed in I Bolton etaL I J2010h . The spectrum is first divided 
through by apower-law F v oc v~ 0,5 , normalised at 1280 (1 + z) A, 
and the Lya emission line is then fitted with a slowly varying 
spline. It is difficult to fit the continuum over the forest at these 
redshifts due to the low flux levels, however in the proximity re- 
gion the transmitted flux maxima will be nearer to the continuum, 
and so is estimated to be within ~ 20 per cent of the correct value 
over the region of interest. 



2.2 Simulated spectra 

Simulations of the IGM at high redshift were used to test the 
method and to explore potential sources of systematic error. The 
simulations are the same as those in iBolton & Beckerl J2009h 
and were performed using a customised version of the paral- 
lel Tree-SPH code GADGET-3, an updated ve rsion of the pub- 
licly available code GADGET-2 dSpringelll2005l) . The simulations 
assume the cosmological parameters (h, Sl m , SIa, fii,/i ,as) = 
(0.72, 0.26, 0.74, 0.024, 0.85), have both dark matter and gas com- 
ponents, and were started at z = 99 with initial conditions gen- 



erated with the transfer function of lEisenstein & Hul dl999l). Each 
simulation uses the UVB model of lHaardt & Madaul J200lh with 
contributions from both galaxies and quasars, and is switched on 
at z — 9 and applied in the optically thin limit. Our fiducial run is 
a 20/2T 1 (comoving) Mpc box with 2 x 400 3 particles (Model C). 
Simulations with other box sizes and particle numbers were used to 
test the dependence of the environmental bias due to an enhanced 
average density on the mass of the host halo, and are summarised 
in Table|2] 

Simulated Lya spectra were constructed from line-of-sight 
density, peculiar velocity, neutral H 1 fraction and temperature 
fields. For the main analysis this was done for 1024 random sight - 
lines drawn parallel to the box boundaries (e.g. lTheuns et aljl998l) . 
with outputs at z = (2,3,4,5,6). Each sightline is 1024 pixels 
long. Proximity zones were introduced into the spectra by modify- 
ing the neutral H 1 fraction (in real space) with the ionising intensity 
falling off as 1/r 2 , before convolving with the other fields to derive 
the optical depth (in velocity space). The spectra were convolved 
with a Gaussian with FWHM equal to the velocity resolution of the 
instrument that was being modelled (6.7 km s -1 for HIRES and 
13.6 kms -1 for MIKE), before being resampled at the instrument 
pixel resolution (2.1 km s _1 and5.0kms -1 respectively). Gaus- 
sian distributed noise, as well as other imperfections, could then be 
added. Additionally, sightlines were drawn through the most mas- 
sive haloes for the analysis described in Section [3.4.2l The haloes 
were identified using a friends-of-friends algorithm with a linking 
length of 0.2. For the analysis in the rest of this paper sightlines of 
40/i _1 comoving Mpc were used. For Models A-C several random 
sightlines were combined for this. 



3 ANALYSIS 

3.1 Proximity effect formalism 

As mentioned in Section Q] the ionising flux from a quasar will 
locally dominate over the UVB in setting the ionisation state of 
the IGM. This leads to increased transmission near to the quasar 
compared to the transmission in the forest (the 'proximity effect'). 
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Figure 1. Top panels: Simulated sightlines at z = 6 (left) and 2 = 5 (right). Each sightline is 40/i — 1 comoving Mpc long, with R cq = 10 proper Mpc and 
S/N = 20. The resolution for the 2 = 6 spectrum is equal to that of HIRES, whilst the spectrum at z = 5 has resolution equal to that of MIKE. The red line 
is the simulated error array. Artifacts such as bad pixels and skyline residuals have also been introduced, hya at the redshift of the quasar is on the far right 
of each panel. Middle panels: The spectra smoothed according to their noise properties (see Section l3~3l , and with optical depths altered such that AF = 
(see Equation llOt . The grey regions mark parts of the spectra that were masked out automatically (primarily skyline residuals). Bottom panels: The original 
Lya forest for these sightlines. The similarity between these spectra and those in the middle panels provides a good check that the method is working correctly. 




8150 8175 8200 8225 8250 8275 7050 7075 7100 7125 7150 
Wavelength (A) Wavelength (A) 



Figure 2. Same as Fig.[T] but for a small sample of the observed spectra. Top panels: The normalised spectrum of SDSS 10836+0054 at z = 5.810 taken 
with the HIRES instrument (left), and the normalised spectrum of SDSS J2225-0014 at z = 4.886 taken with MIKE (right). The length of each spectrum 
corresponds to AOh^ 1 comoving Mpc, and the red line is the observed error array. Lya at the redshift of the quasar is on the far right of each panel. Bottom 
panels: The smoothed spectra with optical depths modified assuming an i? cq of the size quoted in Table|5] The grey regions mark parts of the spectra that 
were masked out automatically. 
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Figure 3. Our adopted fit to the evolution of T e g with r edshift. The blu e 
triangles are from measures of the flux in the forest from Songaila j20O4h . 
whilst the red diam onds are the b i nned v alues from lFan et alj j2006ah . The 
solid line gives the iBecker et al 1 120071) relation for the evolution of T B g 
with redshift, which we adopt for this study. The la error in -r c fj is marked 
with the dashed lines and was taken to be 10 per cent at z = 4, 15 per 
cent at z = 5 and 30 per cent at z = 6, with the error at other redshifts 
calculated by quadratic extrapolation. 



BDO first translated this increased transmission (through a reduc- 
tion in the number of strong absorption troughs) into a measure 
of the photoionisation rate of hydrogen caused by the UVB, Tbkg- 
More recently, measurements have been made of the proximity ef- 
fect using flux transmission statistics rather than line-counting as 
done by BDO. The average transmission measured close to the 
quasar is thereby comp ared to that of the average Lya forest (e.g. 
iLiske & Willigerl200ll) . 

The optical depth, r, is related to the normalised flux, F = 
e~ T . In the simplest model where the quasar lies in a typical region 
of the IGM, and neglecting all motion of the gas and any temper- 
ature gradients that may exist as one approaches the quasar, the 
optical depth at any point can be described as 



Tforest[l + w(r)]" 



(1) 



where Tf orcst is the optical depth that would be measured in the 
absence of the quasar (i.e. the typical value for the forest at that 
redshift), and 



w(r) = 



r q (r) 



(2) 



Here, r q (r) is the Hi photoionisation rate of the quasar at proper 
distance r, and Tbkg is the H i photoionisation rate of the UVB, 
assumed to be spatially uniform for a given redshift. The distance 
from the quasar, r, is approximately 

We further define a characteristic length i? cq to be the distance from 
the quasar where the photoionisation rate from the UVB equals that 
from the quasar, i.e. r q (i? eq ) = Tbkg- With this definition, Equa- 
tion[T]can then be rewritten as 



Tforcst 



1 + 



(4) 



It must be emphasised that R cq is different to proximity re- 
gion sizes as defined in the literature. At high redshift the proximity 
region 'size' is typically defined to be the maximum extent of the 




-0.005 



-0.010 



Figure 4. Top panel: The variation of A_F with log(rbkg) of a simulated 
quasar spectrum at z = 5. The dotted line marks AF = and the dashed 
red line shows the estimated value of log(rbkg). For this particular input 
spectrum A_F = at log(Fbkg) = —12.206. The curve asymptotes as the 
mean flux becomes zero. Bottom panel: Same as the top but for a simulated 
quasar at z = 6. The curve is steeper, and so if the mean flux has been 
overestimated by some factor, this will result in a smaller change in the 
estimated log(r b k^) at z = 6 than at z = 5. 



enhanced transmitted flux (e.g. the first point at which the transmis- 
sion drops to 0. 1 in the spectrum wh en smoothed with a 20 A filter, 
iFan et al. I l2006al : ICarilli et aljkoich . but this is an observational, 
rather than a physical quantity. In this paper we choose to define 
the proximity region size as the scale length out to which the ionis- 
ing flux from the quasar dominates over that from the background, 
i.e. where w = 1, As such any comparison between our proximity 
region sizes (values of R cq ) and those found by other methods for 
the same quasars should be made keeping this difference in mind. 
For further discussio n on the distribution of pr oximity region sizes 
at high redshift see Irvlaselli et ail J2009h and iBolton & Haehneltl 
d2007al) . 

A measurement of the UVB intensity can be expressed as a 
value of i? cq if the flux of ionising photons from the quasar and 
its fall-off with distance are known. We assume a 1 /r 2 fall-off, al- 
though possible deviations from this are discussed in Section [3.4.3| 
For this the flux at the Lyman limit, f VQ , needs to be determined. 
At the redshifts of the quasars discussed in this paper insufficient 
flux is transmitted through the forest for direct measurement of 
fv , and so we extrapolate the continuum flux (measured at a rest 
wavelength of 1280 A as described in Section |2~Tt by assuming a 
power law relation of the form /„ ~ v~ a . The value of a used in 
this paper is 1.61 ± 0.86 and is based on the z < 2 radio-quiet 
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8500 



SDSS J001 1 + 1440, z - 4.967 
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Figure 5. Top: The upper panel shows the normalised spectrum of SDSS J0818+1722 at z = 6.02, whilst the lower panel is a simulated spectrum at z = 6 
with the same R cq and T^kg as presented in Table|5] as well as similar noise properties. The red solid line is the normalised error spectrum, the black dashed 
line is the derived value of -R cq , and the dotted blue line is the expected drop off in the mean normalised flux, hya at the redshift of the quasar is on the far 
right of each plot, and the simulated spectrum has been shifted so that it is at the same redshift as the observed spectrum. Bottom: The same as above but for 
SDSS J001 1+1440 at z = 4.967 and the z = 5 simulation. Again there is a strong similarity between the two spectra. The value of i? cq is smaller at z = 5 
than at z = 6 because the intensity of the UVB is higher, and so the ionising flux from the quasar dominates over the background out to a shorter distance. 



quasar sample of iTelfer et all d2002h . who measured this index in 
the range 500 to 1200 A for 39 individual AG N. The error quoted 
is the RMS scatter of that sample. By contrast. Iscott et all ( 12004) 
found a harder in dex (a = 0.74) in their z < 0.67 sample. The 
IScott et all ( 2004) sample, however, covers a lu minosity range tha t 
is an order of magnitude lower than eith er the Telfer et al] j2002h 
sample, or the quasars analysed here. The lTelfer et al. ( 20021) mean 



value and dispersion should therefore be the most appropriate for 
this study. Using f VQ , the luminosity of the quasar at the Lyman 



limit, L„ , is calculated as 



4ndi 



fr 



(5) 



' (1 + *,) ' 

where is the luminosity distance to the quasar. The quoted errors 
in L„ take into account the error on a, the error in g!l (from the 
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Table 3. The contribution to the error and systematic shift in log(rbkg) from a variety of causes estimated from 1000 simulated spectra at z = 5 and z = 6, 
each with either R oq = 10 or 5 proper Mpc. The errors from specified properties were found to have a nearly Gaussian distribution in log(rbkg). Each 
property, Xi, contributes an error a Xi and a systematic shift, t Xi (both measured in dex). The two instrumental resolutions represent those of MIKE and 
HIRES at z = 5 and 2 = 6 respectively. For the thermal proximity effect the temperature of the gas in the closest 5 proper Mpc has been raised by 10 4 K, 
similar to if the quasar had ionised Hell (see Section [3.4.4) . The total error is calculated by adding the individual errors in quadrature, whilst the total shift 
is the sum of the individual shifts. For comparison, we also simulated spectra with all possible sources of error included simultaneously ('model'). The total 
error, as well as the total systematic shift, are similar to those expected from combining the individual effects. 
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error in the redshift), and the error in the apparent magnitude of the 
continuum, miaso- For quasars with z q < 5.5 the error in mi2so is 
calculated from the error in the measured continuum flux for SDSS 
spectra as described in Section [2"Tl whilst for those with z > 5.5 it 
is taken to be the same as the error on the SDSS z-band photometry. 
For a given distance from the quasar, r, the Lyman limit flux density 
is 



F? (r) = 



47rr 2 



(6) 



The photoionisation rate of H i (in units of s 1 ) by a source of 
UV flux is given by 

AkJ(v)<t H i(v) 



r = 



hv 



-du 



(7) 



where J(v) is the intensity of the source, aniiv) is the ionisation 
cross-section of neutral hydrogen, and h is in this case Planck's 
constant. By definition J bkg (^o) = J^iyo) at r = R aa , and so 
om{v) = 6.3 x 10~ 18 (^o/^) 2 ' 75 cm 2 fcirkman & Tvtleri 



I2008L note that the exponent is often approximated as 3) and inte- 
grating Equation [7] for the photoionisation rate by the background 
gives 



bk g 



9.5 x 10 8 F%(R C 



(8) 



(a + 2.75) 

in units s _1 , where F®(R eq ) is the Lyman limit flux density in 
erg cm -2 evaluated at a distance R cq away from the quasar. There- 
fore, using Equation [6] Tbk g can be expressed as 

9.5 x 10 8 L„ 
bkg - (a + 2.75)47^ ' ( ' 

where L„ is in erg s _1 Hz -1 and R cq is in cm. 
3.2 Measurement method 

The value of Fbkg was inferred for each quasar by increasing the 
optical depths of each pixel by the expected effect of the quasar, 



using Equation [4] until the mean flux across the spectra was the 
same as that expected from the Lya forest at that redshift. Details 
of this procedure are given below. 

The analysis was carried out on the section of the spectrum 
immediately bluewards of the Lya line and 40/i~ 1 comoving Mpc 
(~ 8 — 9 proper Mpc, or Az ~ 0.10 - 0.15) in length. Each 
pixel is converted from a normalised flux into an optical depth us- 
ing t — — \n(F). For a given quasar luminosity and trial value of 
Tbkg, Equation[9]can be inverted to give a trial value of R cq . Using 
this trial value of R eq , the factor [1 + w(r)] -1 can be calculated 
for every pixel in the spectrum. The pixel optical depths are then 
adjusted, before being converted back into fluxes. This process is 
continued until the section of the spectrum resembles the Lya for- 
est at the same redshift. The most robust statistic for determining 
this was found to be the difference between the mean flux in the 
spectrum and the mean flux of the Lya forest, Forest • The differ- 
ence between the two, AF, was defined as 



AF = (F aew ) - F lo 



(10) 



where F new = e~ Tnow and r ncw = t[1 + cj(r)]. The value 
for which AF — we took as our estimate of Fbkg- The 
mean flux, {F nsw }, was calculated as a weighted mean, with each 
pixel weighted by the inverse variance of the flux. For a given 
original error in the transmission, <tf, then <7F ncw = o\f[1 + 
UJ ( r )]~ 1 -Pnow /F. The weighted mean flux was then calculated 

as (F ncw ) = J2 (-Pnew/tr| nsw )/Xi (V°iw)- In order t0 st0 P 
<tf 110W from reaching zero and strongly biasing the weighting, a 
noise 'floor' was used such that crF ncw ^ 0.01. 

Trial values of log(Fbkg) ranged from —14 < log(Fbkg) < 
— 11 and the iteration was done until a precision of 0.001 in 
log(Fbkg) was reached. Fig. Q] shows the effect of the result of 
this iteration on simulated sightlines at 2 = 5 and 2 = 6. The 
'recovered' Lya forest (middle panels) is strikingly similar to the 
actual Lya forest in those sightlines (bottom panels). Two observed 
spectra are presented in a similar fashion in Fig. [2] The bottom pan- 
els show the spectra after smoothing to mitigate the effect of pixel 
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Figure 6. The deviation of the estimated value of log(I\k g ) from the input value for 10 000 simulated spectra at z = 6 (left) and z = 5 (right) with input 
proximity regions of 10 proper Mpc in size. The spectra had a fixed S/N = 20, and the resolution of either the HIRES (top) or MIKE (bottom) spectrographs. 
The dotted black line marks the input value, whilst the blue dashed line marks the mean of the simulated data set. The spread is approximately Gaussian in 
log r with a small systematic offset caused by noise in the spectra. Each offset in the estimated Tbkg corresponds to an overestimate of the input value. The 
width of the Gaussian and size of the systematic offset varied with redshift, S/N, luminosity, and the resolution. Detailed modelling of the effect of each of 
these parameters was carried out to calculate the statistical error and bias in the estimated T^kg for both individual quasars and a grouped sample. A detailed 
eiTor budget is presented in Table[3] and more figures are presented in the Appendix. 



noise (see Section [33t and after the optical depths have been mod- 
ified using the estimated value of R cq presented in Section l4~71 

The assumed value of Ff orC st is calculated at each redshift 
from the fit to me asured optical depths of the Lya forest in 
iBecker et al. I d2007t) . The fit is reproduced along with some ob- 
served values of r e ff from the literature in Fig. [5] It assumes a 
simple evolution of a lognormal distribution of optical depths and 
matches the observed values excellently from 2 < z < 6. Above 
z — 6 there are few measurements of the flux in the forest and so 
it is unknown if this relationship still holds. For the z > 6 quasars 
in this paper it was assumed that it does, although for large values 
of effective optical depth (r > 6), a substantial change in r cor- 
responds to only a very small (absolute) change in Ff orcst . Only 
a small change in R cq , therefore, is required for AF to again be 
zero. This is apparent in Fig. [4] which shows an example of the 
results for one of the simulated spectra at z = 5 and at z — 6. If 
the mean flux is overestimated by a factor of 5, for example (corre- 
sponding to the dotted line being down at —0.10 and —0.010 in the 
top and bottom panels respectively) then the shift in the estimated 
logfTbkg) is much greater at z — 5 (~ 0.6 dex) than at z — 6 
(~ 0.3 dex). The results from the z > 6 quasars should therefore 
be rather robust to significant uncertainties in the mean flux. 

Probable skylines are identified as regions with of > 2 (<7f), 
where the mean is defined over the whole spectrum. Those pixels, 



plus 5 pixels either side, are masked out. The data is also smoothed 
with a boxcar with a smoothing window that is proportional to 
the amount of noise in the data (window width = 10[(ctf) /0.05] 
pixels), in order to balance between a smoothly varying function 
whilst maintaining some of the contrast and resolution in the prox- 
imity region. Regions with a smoothed flux below zero were treated 
as though they were positive, but remembering the sign, such that 
Ficw = -e~ ( ~ ln|F|)[1+ '* J(r), . Smoothed fluxes that were greater 
than or equal to 1 (i.e. at the continuum) were set to be equal 
to 0.99. These effects, designed to limit the effects of observa- 
tional artifacts, introduce a bias such that as the noise increases 
the method will systematically overestimate the UVB. These bi- 
ases were therefore extensively modelled using the simulated spec- 
tra (see Appendix). 

There are several advantages of our method over other recently 
presented methods to measure the intensity of the UVB from the 
proximity effect. One of the key attributes of our method is that 
the intrinsic distribution of Lya optical depths need not be known 
a priori. This means that we are not directly dependent on numer- 
ical simulations. Another key advantage is that a slope, j3, of the 
column density dist ribution of neutral hydro gen is not assumed, as 
was required in e.g. iDall' Aglio etai] d2008h . We believe that with 
our method the largest remaining uncertainty is the evolution of T e g 
at z > 6. 
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Figure 7. Demonstration of the impact of Lyman limit systems (see Section l3.4.3l . Top: A standard simulated spectrum at z = 6 with i? cq = 10 proper Mpc 
and S/N = 20. The thick black dashed line marks this input R cq whilst the dotted red line shows the one estimated by our method. The error bars correspond 
to the statistical error in the estimated R cq , roughly 15 per cent. The thin black dashed line is the expected drop in normalised flux for the input R C q- The 
method has, within error, recovered the input R cq . Middle: The same spectrum as above, but with a Lyman limit system with corresponding Lya absorption 
at the location of the blue arrow. The region of enhanced transmitted flux is slightly truncated, and so the estimated R cq is slightly smaller than in the case of 
no Lyman limit system, corresponding to a higher estimated value of Tbkg- This is a well known bias in proximity effect measurements. Bottom: The same 
spectrum as in the top panel, but with the Lyman limit system much closer to the quasar, causing substantial shortening of the apparent size of the proximity 
region. There is again an associated decrease in the estimated H e q, and thus a higher estimated r bkg . However, the estimated values of R oq are substantially 
more robust than would be expected from the extreme shortening of the apparent extent of the enhanced transmitted flux. 



3.3 Statistical accuracy of the method 

In order to test the accuracy of the method, simulated spectra were 
created using random sightlines through the simulations. The test 
spectra were all created using Model C at z — 5 and 6, and were 
40/i -1 (comoving) Mpc long, constructed by joining together mul- 
tiple sightlines. In our standard mock spectra the proximity region 
was assumed to have a size of 10 proper Mpc. This is in gen- 
eral agreement with the deriv ed values from t he dat a (see Sec- 
tion [4J3 and that estimat ed bv lWvithe & Loebl d2004) . Assuming 
a lHaardt & Madaul fcOOlh UVB this corresponds to a quasar with 
M1450 = -27.10 (AB) at z = 6 and M 145 o = -27.47 (AB) at 
z — 5. The noise was assumed to be Gaussian with a signal-to- 
noise (S/N) of 20, close to the average value of the S/N in the 
data. 

In order to better mimic the real data, a variety of artifacts 
were introduced into the simulated spectra. Sky line residuals were 
added with an associated peak in the error spectrum, as well as bad 
pixels for which there was no associated effect on the error spec- 
trum. Overall the simulated spectra are qualitatively very similar to 
the real data, as demonstrated in Fig. [5] 

We applied our analysis to artificial spectra with similar S/N 
to the data. On average the method estimated the correct value of 



Tbkg, with the distribution of values approximately Gaussian in 
log(rbkg). Several factors were thoroughly investigated for their 
effect on the size of the errors, and any bias on the estimated UVB 
intensity. We investigated the effects of peculiar velocities in the 
gas, Vp C c, which will introduce distortions in redshift space to the 
expected radial profile of the transmitted flux. This was done by 
re-running our analysis on spectra generated with peculiar veloci- 
ties set to zero. Similarly we examined the impact of sightline-to- 
sightline variations in the mean flux away from the global mean 
flux due to local fluctuations in the density field. For this we ran 
the analysis on spectra from which the mean flux of that sight- 
line was known (equivalent to if the quasar was removed perfectly), 
and compared it to when simply assuming a global mean flux. The 
spread was also dependent on the noise levels in the simulated spec- 
tra, the redshift, the instrumental resolution, the luminosity, and 
the error in rarest (which varied with redshift). This analysis is 
shown in more detail in the Appendix. Each contributor to the er- 
ror (see Table |3) was varied between reasonable limits for 1000 
simulated spectra at both z = 5 and z = 6, and linearly extrap- 
olated to other redshifts. This meant that the contribution to the 
error budget and associated bias from all the variables could be cal- 
culated for the parameter space and redshift range present in the 
data. Using this analysis an overall statistical error and bias were 
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determined for each individual quasar. For demonstration, the de- 
pendence on S/N and instrumental resolution is shown in Fig. [6] 
As expected, if the S/N decreases then the scatter in log(rbkg) in- 
creases. There is also a systematic bias, due to the smoothing of the 
spectra. Averaging over fluxes is not equivalent to averaging over 
optical depths, due to the non-linear relation between the two, and 
so the more the spectrum is smoothed the more the optical depths 
are underestimated. This effect is strongest for low fluxes (high op- 
tical depths). Consequently, smoothing (particularly at high red- 
shift) increases the transmission in the forest so the fall-off in mean 
flux is not as steep. This corresponds to a smaller R cq , and so ulti- 
mately log(rbkg) is overestimated. Since MIKE has a lower instru- 
mental resolution than HIRES there is an additional bias brought in 
due to a similar effect. Indeed any effect that smoothes the spectra 
by averaging over pixel fluxes will bring in a bias of this nature. 
Table [3] breaks down the sources of scatter and systematic biases 
based on 1000 sightlines that have similar values to those in the 
data at z = 5 and z = 6. At both redshifts the errors are dominated 
by the error in r e ff, and the errors in redshift are far more important 
than the errors in the luminosity. Even though the error in r e fs is 
symmetric (it is approximated as a Gaussian), the resulting bias is 
asymmetric, and becomes more so at high t c h (i.e. high redshifts), 
as the equivalent distribution of Ff orca t also becomes highly asym- 
metric. 

As we noted in Section [37X1 we assume a spatially uniform 
UVB. However, shortly after reionization has completed the mean 
free path of ionising photons is short enoug h that significan t 
spatial variations in the UV B may exist (e.g iLidz et al] |2007|) . 
iMesinger & Furlanettd d2009h show however that variations in the 
density field at z ~ 5 — 6 dominate over spatial variation in the 
UVB, and that even just after reionization assuming a uniform UVB 
underestimates Fbkg by at most a few percent. We thus do not try 
to correct for this uncertain but small effect. 

We also performed a joint analysis of several sightlines simul- 
taneously. These were modelled in a similar way, with the lumi- 
nosities, noise properties and instrumental resolution equal to that 
of the component spectra in each bin being represented. From this 
a statistical error and bias could be determined for each bin. 

The effect of other sources of systematic bias, such as Lyman 
limit systems (LLS, see Section [3.4.3b . were also considered. These 
LLS are absorption features that are optically thick (r > 1) to Ly- 
man limit photons, and so prevent the quasar from ionising as large 
a volume of the IGM as in their absence. Consequently they can 
truncate the extent of enhanced transmission in a spectrum. Fig. [7] 
demonstrates the case where a Lyman limit system has shortened 
the apparent proximity region. As the shortening becomes more 
severe, the method underestimates the true proximity region size 
(and thus overestimate the value of log (Fbkg)) as the LLS modifies 
the assumed 1/r 2 fall-off from the quasar. However, our method is 
more resilient to this effect than previous proximity effect determi- 
nations of the UVB in the literature. Consequently it will produce 
on the whole more accurate values of UVB, less susceptible to the 
effects of Lyman limit systems. This is discussed in more detail in 
the following section. 



3.4 Details of potential systematic biases 

Whilst the pr oximity effect can be used on both individual 
sightl ines (e.g. Iwilliger et alii 1994 ICristiani et~ai]| 19951 : iLu et all 



UVB estimates from the proximity effect can be considered robust. 
The proximity effect is essentially measuring the direct signature 
of an enhanced ionisation field near the quasar by noting a reduced 
fraction of H 1 extending out across a region of a few physical Mpc. 
Therefore any conversion from this physical size to a value for the 
UVB will be highly sensitive to any assumptions made about the 
quasar and its environment. There are three main systematic uncer- 
tainties in this conversion. 

First, it is assumed that the quasar turned on a sufficiently long 
time ago that its proximity region is in photoionisation equilibrium 
with the IGM, and that it has been at its current luminosity over 
a similar timescale. The unknown level of variability in quasar lu- 
minosity on the order of 10 4 years (the timescales required for a 
region of the highly ionised IGM to reach ionisation equilibrium) 
means that the measured proximity region size may not correspond 
to a region in ionisatio n equilibrium, and so any measurement of the 
UVB may be biased Jpentericci et al J 1 20021 : Ischirber et ai]|2004l : 
iDall' Aglio et alj2008h . 

The second problem is that in order to calculate an accu- 
rate size of the proximity region, an accurate systemic redshift 
for the quasar is needed. Redshifts determined from broad high- 
ionisation lines are underestimates of the systemic redshift (e.g. 
iRichards et alj|2002h . and an underestimate in the systemi c red- 
shift will cause the UVB to be overestimated tespevill993l) . It is 
worth noting, however, that an accurate systemic redshift can be 
determined from low ionisation lines, and so for some spectra this 
effect can be avoided. Also it should be noted that at high redshift 
the extent of enhanced transmitted flux will be larger than even the 
largest error in systemic redshift present in the data. 

Finally, in the standard proximity effect analysis by BDO, 
the IGM within the proximity region is presumed to have the 
same density distribution as the IGM outside it. Quasars are how- 
ever hosted by massive galaxies which are expected to reside 
in an environment whi ch has higher than a verage density out to 
rather large distances jGranato et al .1 120041 : iFontanot et al.l 1 20061 : 
Ida Angela e7al]|2008l) . In that case, the proximity region would be 
smaller than expected for a region of the same size with an av- 
erage density close to the global mean densit y (jD'Odorico et al 
2008 ), and thus the UVB may be overestimated JLoeb & Eisenstein 



1996) and i n a statistical sense acr oss many spectra (e.g. BDO; 



Cooke et alii 19971 ; IScott et alj|200d) to help constrain the redshift 



evolution of the UVB, several issues must be addressed before the 



1995). This environmental bias should be worst for the most lu- 
minous quasars which presumably lie in some of the most ex- 
treme overdensities at a given redshift (e.g . IPascarelle et alj|200ll : 
lAdelberger e"t~ai1l2003l ; iKim & CroftH2007l) . As already discussed 
this is believed to be the main reason for the discrepancies be- 
tween estimates of t he UVB using the pro x imity effect and the flux 
decrement method 1 Rollinde et al. 1 120051 ; iGuimaraes etafl 120071 : 
iFaucher-Giguere et alj2008l) at redshift z ~ 2 - 4. 

Another environmental effect that is often neglected is the 
consequence of the quasar heating the surrounding IGM via ion- 
isation of He 11 . lBoltonetai] l l2010l) found that the gas within ~ 5 
proper Mpc of the z « 6 quasar SDSS J0818+1722 was ~ 10 4 K 
hotter than that presumed for the general IGM at that redshift. 
Higher temperatures within the proximity region will lead to more 
transmission, and so the UVB might then be underestimated. 

All of these effects were investigated with the simulations at 
z = 5 and z — 6 to try and quantify how important they might be 
in the real data, and so each of these points will now be discussed 
in detail. 
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Figure 8. Mean density profiles around the most massive haloes in our 
simulation boxes. Top: The solid line is the average density field in the 
sightlines that start in haloes of masses greater than 10 12 Mq in Model F 
at 2 = 6, plotted in proper Mpc. The dotted lines mark the ±1ct deviations, 
and the mean density in the simulation is marked with the horizontal dashed 
line. Any significant overdensity is restricted to a region of ~ 1 Mpc in 
these haloes, and the overdense region will be even smaller for lower mass 
haloes. Bottom: The same as above but for haloes with mass greater than 
10 13 Mq at z = 3. These are essentially the same haloes as those at z = 6, 
but have become more massive through hierarchical growth. The overdense 
region is larger, extending out to ~ 3 Mpc, with evidence for smaller haloes 
clustered nearby (the small peak at ~ 5 Mpc). If the UVB is higher at this 
redshift then proximity region sizes will decrease, and may be comparable 
to the extent of the overdensity. In that case measurements of the size of the 
proximity region could depend strongly on the mass of the quasar host halo. 



3.4.1 Effect of luminosity and redshift errors 

Investigation with both the real and simulated spectra showed that 
the measured error on the UVB is dominated by errors in the red- 
shift and the effective optical depth of the forest, and not by errors 
in the quasar luminosity. Even random errors in luminosity of up to 
40 per cent had only a comparable effect to small errors (~ 0.01) in 
redshift. Of course, if the variability in quasar luminosity is not sim- 
ply an extra source of random scatter, and in fact is instead system- 
atic (for example, if the quasars are only observed in their brightest 
stages) then there will be a systematic bias as well as an additional 
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Figure 9. Top row: The ratio of the estimated l\kg (corrected for system- 
atic shifts due to noise) to the input T^kg for sightlines starting in the 500 
most massive haloes in the fiducial simulation, Model C. The proximity re- 
gion in each spectra is 10 proper Mpc. On the left results are shown for 
2 = 6 and on the right z = 5. The black dashed line is for zero bias, whilst 
the red line marks the true average value estimated. For these haloes there 
is a systematic overestimation of < 0.1 dex, but apparently no dependence 
on halo mass. Middle row: Same as above, but for Models B (black) and E 
(blue). Model B has been resampled to the same pixel size as Model E. Bot- 
tom panel: Same as above but for Models A (black), D (blue) and F (green). 
Due to the low resolution and small box size in Model A only 100 haloes 
could be identified at z = 6 and 250 haloes at z = 5. Models A and D have 
been resampled to the same pixel size as Model F. Lowering the mass res- 
olution introduces a systematic overestimation in log(rbkg). This is due 
to the fact that the voids are under-resolved and so the proximity regions 
appear systematically smaller. The larger pixel sizes used in the lower mass 
resolution simulations also smooth the transmitted flux, so the proximity 
region is again underestimated (see Section [33t . This effect is the stronger 
of the two. With these offsets taken into account, there is still no significant 
increase in l\k g with halo mass, even in the bottom panel, which includes 
an 80/i -1 comoving Mpc box with haloes of mass > 10 12 Mq. 



random error. The method ultimately is sensitive to the gradient of 
the fall-off in the mean transmitted flux throughout the spectrum, 
and so will give the same estimate for R cq regardless of the quasar 
luminosity. Thus, if all the quasars are radiating at a systematically 
brighter luminosity than the one that established their proximity re- 
gion, then the values of Tbkg will have been overestimated by the 
same amount, since they are linearly related in Equation[9] 

Several of the quasars in this paper have a very accurate sys- 
temic redshift from either CO or Mg n emission lines dCarilli et al.l 
l2010l) . and the high resolution spectra themselves allow for a rea- 
sonably accurate determination of the onset of the Lya forest. Con- 
sequently almost all the redshift errors are ^ 0.01, corresponding 
to a size much smaller than the proximity region. We therefore be- 
lieve that these two possible sources of error have been suitably 
dealt with. 
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3.4.2 Effect of quasars lying in overdensities 

One of the more contentious assumptions used in proximity re- 
gion measurements is that the environment of the quasar within the 
proximity region is similar to that of the general IGM. The quasar 
should lie in an overdense environment, and therefore could strug- 
gle to ionise as large a volume as it would in the general IGM, thus 
reducing the apparent size of the proximity region. As discussed al- 
ready there have recently been a number of studies suggesting that 
proximity effect mea surements may indeed suffer from a substan- 
tial bias due to this (Rolli nde et al, I l2005l : iGuimaraes etai]|2007l ; 
iFaucher-Giguere et alj2008h . 

We have therefore tried to estimate how much such a bias may 
affect our measurements, which are at significantly higher redshift 
than previous studies. As a first attempt we looked into the effect 
of placing the quasar at the centre of massive dark matter haloes in 
our fiducial simulation and simulations with different box size and 
mass resolution. 

For this we have chosen sightlines through the largest (by to- 
tal mass) 500 haloes in Models B-F in 3 perpendicular directions 
through the box. Due to the low resolution and small box size of 
Model A, only 100 haloes could be identified at z — 6 and 250 
at z = 5. For each halo the sightline was adjusted so that it be- 
gan at the centre of a massive halo. Since the simulation box is 
cyclic the sightlines needed to be split in half so the sightline only 
passed through the overdensity once, and consequently for each 
halo 6 sightlines could be drawn (i.e. leaving the halo from both 
sides in 3 perpendicular directions). The halo sightlines were then 
spliced with random sightlines so that the total length of each spec- 
trum was 40/i _1 comoving Mpc. The proximity region size was 
chosen to be 10 proper Mpc, as in previous tests with the simula- 
tions. The results for the haloes in the fiducial simulation, Model C, 
are shown in the top panel of Fig. [9] There is no significant change 
as the mass increases, with a Gaussian spread in log(Tbkg) in each 
bin of similar width to that from the random sightlines (Fig.[6jl. 

The box size for Model C is, however, only 20ft _1 comov- 
ing Mpc, which limits the maximum halo mass to just over 4 x 

10 11 M@. Luminous quasars are likely to reside in host hales more 
massive than this. The best observational constraints are inferred 
from the clustering analysis of quasars combined with predictions 
from ACDM simulatio ns and suggest host halo masses in the range 

10 12 - 10 13 M© (e.g. iBonoli et al]|2009h . It is important to note 
here that the very luminous quasars often used for proximity effect 
measurements are more luminous than those available for cluster- 
ing analysis. It is thus perhaps not surprisin g th at the proximity 
effect studies studies bv lRollinde et al .1 d2005l) and lGuimaraes et ahl 
d2007l) suggest that at 2 ~ 2.5 — 4 the host haloes of these very 
luminous quasars are even more massive. Fortunately, at the red- 
shifts we consider here haloes more massive than 10 13 Mq are not 
yet expected to have formed in large enough numbers to be cred- 
ible candidates to host even the most lu minous observed quasars 
dSpringel et all 120051 : ISiiacki et al.ll2009h . In order to probe over- 
densities around haloes in this mass range larger boxes are needed 
than we have considered so far, and so simulations with 40ft _1 and 
80h~ comoving Mpc box size were also investigated. These were 
then compared with 20h~ 1 Mpc boxes with equal mass resolution. 
Model B was compared to Model E (see middle panel in Fig. [9]l, 
and Model A to Model D and Model F (bottom panel). Sightlines 
through Models A, B and D were resampled so that the pixel size 
(and mean forest flux) was the same as the largest simulation to 
which they were being compared. The highest mass bin in Model F 
(> 10 12 M ) still shows no more significant bias (in comparison to 



less massive haloes in Models A and D), seeming to confirm the re- 
sults from Model C that, at least at z = 5 — 6, the mass dependence 
of the bias is rather weak. 

The largest halo mass in the z = 6 simulation was 1O 122 M , 
which may still be nearly an order of magnitude less massive than 
the haloes the luminous quasars used in our study may reside in. To 
get a rough estimate of the effect of more massive haloes we have 
taken the density profile of haloes with M > 10 13 M at z — 3 in 
our simulation (the bottom panel in Fig.(8]( and modified the optical 
depths of spectra created from random sightlines in the z = 6 sim- 
ulation by using r oc A 2 ' 7-1 ), where 7 is from the temperature- 
density relation, T(A) oc A 7-1 , and is taken to be 1.3. The extra 
systematic bias from such a massive halo is shown in Tablef3] The 
effect is small at z = 6, and even though it is more important at 
z = 5, it is still not the main source of bias or main contributor to 
the error budget. Obviously this is just an approximation of how a 
halo this massive would affect the results, as the density profile was 
applied to the spectra artificially and will not be self-consistent with 
the other simulation outputs along the line-of-sight. However, since 
the effect is relatively minor we would not expect it to be a domi- 
nant source of error even if it were modelled self-consistently using 
a halo found in a much bigger simulation box. Since the quasars at 
z — 6 will probably lie in haloes with M < 10 M our estimate 
should be an upper limit of the effect of overdensities, and thus we 
conclude that at z — 5 — 6 the expected overdensity around our 
quasars is not critically important. 

The most likely explanation for the lack of dependence of the 
bias of the UVB estimate on the mass of the host halo is that the 
proximity region is large in comparison to the size of the overden- 
sity. This is demonstrated in Fig. [8] where the overdensity profile 
is shown at z — 6 and z — 3 for the largest haloes in Model F 
(M > 10 12 M and M > 10 13 M respectively). The gas will 
thereby tend to be infalling and the corresponding peculiar veloci- 
ties will further de crease the apparent size of th e proximity region 
in velocity space dFaucher-Giguere et alj|2008l) . In our simulated 
spectra the proximity regions were 10 proper Mpc, which is nearly 
an order of magnitude larger than the overdensity expected in the 
largest haloes at z — 6, and so the properties of the IGM within the 
proximity region should not be strongly biased due to the overden- 
sity. 

As a final point, we note that some clustering of galaxies and 
faint Active Galactic Nuclei (AG N) is expected around luminous 
quasars (e.g. lUtsumi et alj|201fj|) . The intensity of the UVB in the 
vicinity of luminous quasars could therefore be enhanced due to 
other nearby sources. Our knowledge about the clustering of faint 
AGN and galaxies around bright high-redshift quasars is rather 
sparse. We can nevertheless use the magnitu des and posit i ons of 
the Lyman Break Galaxies (LBGs) found bv lUtsumi et al .1 d2010t) 
for a rough estimate. In a field containing a quasar at z = 6.43 
they found 7 LBGs at z > 6.4, based on their colours. Let us 
then calculate L VQ for each LBG from their zr band magnitude, 
assuming they radiate isotropically, are are all at the same red- 
shift as the quasar, and have a ratio of fa r-UV to extreme-UV flux, 
/1500//900 = 22 dShaplev et alj|200oi as measured for z ~ 3 
LBGs. With projected distances of 2 — 5 Mpc the LBGs would all 
lie within the typical proximity zones of our sample if they were 
indeed at the same redshift. Even if the LBGs would have a Lyman 
continuum as blue as we assume for the quasars (i.e. a = 1.61) 
they would increase the size of the proximity region by as little as 
~ 0.1 Mpc, and the UVB would typically be underestimated by 
0.01 dex (3 per cent). There could obviously be a substantial con- 
tribution from more numerous fainter objects below the detection 
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Table 4. Summary of the parameters of the different power laws used to investigate the effect of LLS on both individual and grouped spectra. The number 
of expected LLS in a random sightline through Model C is calculated using Equation 1111 with z max equal to the simulation redshift and z ma x — Zmin 
corresponds to 20h~ 1 comoving Mpc. SC10 parameterise their power law with jVrj evaluated at z = 3.5 such that instead of N(z) = 7Vo(l + z) 7 they have 
N(z) = ./V3.5 [(1 + z)/4.5] 7 . The shifts, e, quoted are the average overestimation in log(rt>k g ) due to the presence of LLS. For the individual spectra, those 
at z = 6 have HIRES resolution, and those at z = 5 have MIKE resolution. For the grouped spectra the shifts correspond to the z ~ 5 and z ~ 6 redshift 
bins, as described in Section l4~Tl The final row represents the results based upon a power law consistent with the 1 u upper limit of the SC10 parameters. 





Redshift range 


N 


^3.5 


7 


"lis 


z = 5 
e (ind) 


e (bin) 


"lis 


2 = 6 

e (ind) 


e (bin) 


SL94 
P03 
SC10 
SC10 la 


0.40 < z < 4.69 
2.40 < z < 4.93 

< z < 6 


n 97 +0.20 

07+ - 13 

u - ul -0.04 

0.15 a 


2.80 ±0.33 


1.55 ±0.45 
9 4^+0.75 

Z ' 40 -0.65 
1 G4+0.36 


0.2169 
0.2810 
0.2441 
0.3022 


0.15 
0.20 
0.17 
0.21 


0.09 
0.12 
0.10 
0.12 


0.3459 
0.5146 
0.4134 
0.5409 


0.17 
0.26 
0.22 
0.27 


0.08 
0.11 
0.09 
0.12 



a Derived from their -/V3.5 and 7 values. 



limit. The faint end slope of the luminosity function of LBGs at 
z ~ 6.4 is not well co nstrained, but if we assume a slope of —1.73 
jBouwens et alj|2010l) down to zero luminosity in the volume of 
the proximity region the UVB is still only underestimated by 0.04 
dex (11 per cent). This is, however, likely to be an upper limit, as 
in reality the 7 LBGs do not have precise redshifts and may not be 
within the proximity region of the quasar. This should outweigh the 
possibility that the galaxies may have a larger value of /1500 / /900 
than we have assumed here. We therefore conclude that the effect 
of nearby LBGs could be noticeable but appears likely to be small. 
Considering the large uncertainties in our rough estimate we did 
not try to correct for this. 



3.4.3 Effect of Lyman limit systems 

Lyman limit systems (LLS) are due to regions of neutral gas that 
are optically thick (r > 1) to Lyman limit photons (A lcs t = 
912 A), which corresponds to a hydrogen column density of 
JV(Hi ) > 1.6 x 10 17 atom cm" 2 . Should a LLS lie in the ob- 
served line-of-sight of the quasar then the apparent size of the prox- 
imity region will be shortened, as th eir abundance limits the mean 

Storrie-Lombardi et all 1 1994 



covering the range 0.40 < z < 4.69, and suggest that for a sight- 
line covering the redshift range z mBX to z m i n that 



free path of ionising photons (e.g 



Miralda-Escude 2003; [Peroux et al.l 



20091) . 



20031 ; iFurlanetto & Mesinged 



Studies counting the number of LLS in spectra struggle at high 
redshift as, due to the ever increasing blackness of the forest, fea- 
tures with r > 1 are difficult to detect. Consequently there are very 
few studies that provide measured LLS frequency at the redshifts 
covered in this paper. LLS can lead to dramatic shortening of the 
region of enhanced transmission (see Fig.fTJl and so if abundant at 
high redshift then they could be a substantial source of systematic 
errors. 

We inserted LLS into o ur simulations following th e method 
presented in Appendix D of iBolton & Haehneltl J2007al) . A den- 
sity threshold, p t hresh, was chosen such that the average num- 
ber of regions in a spectrum with p ^ /Jthresh was the same as 
a given number of LLS expected per sightline. The neutral frac- 
tion within those regions was then set to unity (i.e. they become 
self-shielded). As such, they absorb all the ionising flux of the 
quasar so further out in the spectrum the transmission in the for- 
est is from the UVB alone. The expected number of LLS in a 
particular sightline, wn B , was calculated from the power laws pre- 
sented by Storrie-Lombar di et al. I fl994h .l Peroux et al. IHH), and 
ISongaila & Cowiei fcOld) . (hereafter SL94, P03 and SCI respec- 
tively), where N(z) = dN/dz = N (l + z) 1 . SL94 used quasars 



nils = 



2Vb(l- 



(11) 



where No — 0.27 and 7 = 1.55. For a particular random sightline 
in Model C z max was taken to be the redshift of the simulation and 
z mal — z m i n corresponds to the simulation box size. This meant 
that mis = (0.2169,0.3459) at redshift z = (5,6) respectively. 
From this (using only the random sightlines) the threshold den- 
sities were derived as log(A t h rC sh) = (1.6310, 1.4080) (where 
Athresh = Pthresh/pav) and those regions were then presumed to 
be self-shielded. Similar calculations were also done for the power 
laws of P03 and SC10, and a worst case scenario was investigated 
using the la upper limits of SC10. The expected overestimations 
in log(rbkg) for each of these power laws are presented in Table|4] 
for both individual spectra, and grouped spectra. Consequently, as- 
suming the predictions of SC10 to be the most reliable as this is the 
only study probing the redshift range we are interested in, the effect 
would be to shift the results of the individual quasars down by 0. 17 
(0.22) dex at z = 5 (6), and possibly up to 0.21 (0.27) dex. 

SC10 is the only study into the actual spatial frequency of LLS 
at z = 6. In both SL94 and P03 they found that towards higher red- 
shifts, N(z) starts to evolve rapidly, and so extrapol ations become 
much more uncertain. Indeed lProchaska et al .1 J201LT) note that both 
of those papers are subject to biases that they claim have not been 
adequately compensated for, and so they may have overestimated 
N(z) by up to a factor of 3. Clearly if nn a is a factor of 3 smaller 
then the effect on Tbk g would be greatly reduced. SC10 themselves 
also warn against extrapolating to z > 6 due to the potential rapid 
change in the mean free path of ionising photons if reionization 
is being probed. Taking all these factors into account, and treating 
SC10 as the most reliable power law, we cautiously estimate using 
their la upper limits on nn s that the binned data points in Fig. 1131 
are possibly overestimated by ^ 0.12 dex. 



3.4.4 Effect of quasars heating the IGM in the proximity region 

We have implicitly assumed that the temperature of the quasar 
proximity region is comparable to that of the general IGM. This 
may be not true, particularly prior to He 11 reionization, when the 
ionisation of Hen by the quasar may heat the local IGM (e.g 
lBoltonetal1l2010h . This temperature gradient will cause higher 
mean flux in the proximity region, and so the proximity region 
size will be overestimated and the UVB underestimated (i.e. op- 
posite to all the other environmental biases). We constructed a toy 
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Figure 10. Left panel: The bias corrected estimates of log(Fbkg) for the 15 individual quasars studied in this paper. The errors are calculated as the standard 
deviation in log(Fbkg) from 10 000 Monte Carlo realisations (varying the systemic redshift, luminosity and forest r c g) added in quadrature to the statistical 
error, which was individually calculated for each quasar sightline. Biases similar to those listed in Table[3]have also been removed. A smooth decline in r bkg 
with redshift appears over this redshift range, with a formal correlation coefficient of —0.87. Right panel: The inferred values of R cq for each quasar. The 
error bars were calculated in the same way as those for r bkg . There is a strong increase towards higher redshifts, driven by the declining intensity of the UVB. 
A fit to this increase is shown by th e solid line, and the forma l correlation coefficient is 0.86. The dashed line marks the evolution in the mean free path of 
ionising photons from the formula in ISongaila & Cowidfepioll The size of the proximity region is larger than the mean free path at z > 5.5. The implications 
of this are discussed in Section l4~2l 
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Figure 11. Plot of the deviance from the linear fit to log(r b i cg ) in the 
left panel of Fig. [10] against the quasar luminosity. The dashed line marks 
zero deviation. All of the points are consistent with this line, which suggests 
that the measured values of log(Fbkg) have been adequately corrected for 
systematic biases from quasar luminosities. 

model in the simulations with the change in temperature as a step 
function, such that the gas within 5 proper Mpc of the quasar was 
10 K hotter than that in th e general IGM. This is similar to the re- 
sult of lBolton et all feoiot). assuming gene ral IGM temperatures at 
z = 5 - 6 of ~ 10 4 K teeckeretalfeoiOt) . In individual sightlines 
this caused log(Fbkg) to be underestimated by 0.10 (0.11) dex at 
z — 5 (6), and in the binned data by 0.14 (0.21). Consequently, 
we find that the hotter temperatures of gas close to the quasar could 
cause underestimation in log(rb kg ) comparable to the overestima- 
tions from other environmental effects at the redshifts of our study. 



4 RESULTS AND DISCUSSION 
4.1 Results 

Table [5] summarises the main results of this paper for the indi- 
vidual quasar sightlines. The quoted errors in Tb k g consist of the 
measured error from 10 000 Monte Carlo realisations of the spec- 



tra, varying the luminosity, systemic redshift and forest T e g, added 
in quadrature to the expected sightline-to-sightline scatter as de- 
scribed in Section [33] The spectra with the largest errors in L^ g 
are those with continuum magnitudes measured directly from the 
SDSS spectra. Even though the errors in the luminosity could be 
rather large (~ 30 per cent), the measured error was strongly dom- 
inated by errors in the redshift and r e s- The statistical power of a 
measurement from a single spectrum is rather limited and the er- 
ror on Fbkg is substantial. Fig. [10] shows the estimated log(Fbkg) 
and R cq for our sample as a function of redshift. The data points 
in Fig.[l0]are consistent with a linear fit within the errors. There is 
therefore little we can say about spatial fluctuations in UVB inten- 
sity other than that they appear to be smaller than our m easurement 
errors . We note that if we had used the a presented in lScott et al.l 
d2004l) then the log(Fbkg) values would have been on average 0.22 
dex higher, and the error bars 10 per cent larger. We also note that 
our redshift evolution of R cq is opposite to the redshift evolu tion of 
the proximity region sizes as presented bv lFan et al .1 d2006al) . Note 
that this is entirely due to the different definition of the two sizes 
conside red. For the proximity region size definition from lFan et al.l 
d2006d) we find similar results. 

In Fig.[TT]we subtract out the linear evolution of the average 
log(Fbkg). There is no systematic trend with the luminosity of the 
quasar. This confirms that the systematic shifts induced by the lu- 
minosity of the quasar (see Tableland Fig. lA5l > have been suitably 
corrected for. Fig. [T2] shows a sample of the observed spectra with 
their derived R cq and expected average flux fall-off, in their rest 
frame, as well as denoting the section of each spectrum used for the 
measurement. Judging by visual inspection, none of our observed 
spectra appears to be influenced by LLSs to the degree presented in 
the bottom panel of Fig. [7] It is, however, difficult to determine if 
any of the sightlines contained LLSs far from the quasar redshift, 
as in the middle panel of Fig.[7] SDSS Jl 148+5251 has a very short 
region of enhanced transmitted flux (Fig. 112b and is a possible can- 
didate for being affected by a LLS. 

Our most robust results come from simultaneously fitting mul- 
tiple sightlines. The data were grouped into low redshift (z < 5.5) 
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Table 5. Tabulated r esults for each of the investigated quasar sightlines. Those with a z of 0.002 or 0.007 have systemic redshift taken from the CO or 
Mg II emission lines ICarilli et alJl201oT) while the others use either the Lya + Nv emission line redshift, or (in most cases) use the redshift at which the 
Lya forest appears to begin. Due to the high-resol ution of the spec t ra, errors for these a re 0.01. The absolute magnitudes at rest-frame 1450 A were calculated 
from either the published continuum magnitudes iFan et al. 2001, 2003, 2004, 2006b) or the measured fluxes at rest-frame 1280 A, extrapolated to 1450 A by 
assuming a power law continuum of the form f v oc v~ °' s . 



Name 











10 31 ergs -1 Hz -1 


(AB) 


SDSSJ1 148+5251 


6.4189 


0.002 


CO 


3.154 ±0.935 


-27.81 


SDSS J1030+0524 


6.308 


0.007 


Mgn 


1.721 ± 0.516 


-27.15 


SDSS J1623+31 12 


6.247 


0.007 


Mgn 


1.100 ±0.329 


-26.67 


SDSS J0818+1722 


6.02 


0.01 


Lya forest 


2.156 ±0.649 


-27.40 


SDSS J1306+0356 


6.016 


0.007 


Mgn 


1.792 ± 0.536 


-27.20 


SDSS J0002+2550 


5.82 


0.01 


Lya forest 


2.754 ±0.829 


-27.67 


SDSS J0836+0054 


5.810 


0.007 


Mgn 


3.333 ± 0.998 


-27.87 


SDSS J023 1-0728 


5.41 


0.01 


Lya forest 


1.323 ±0.738 


-26.87 


SDSS J1659+2709 


5.33 


0.01 


Lya forest 


2.516 ± 1.079 


-27.57 


SDSS J0915+4924 


5.20 


0.02 


Lya + N V 


1.744 ±0.830 


-27.17 


SDSS J1204-0021 


5.09 


0.01 


Lya forest 


1.900 ±0.851 


-27.26 


SDSS J001 1+1440 


4.967 


0.005 


Lya forest 


4.518 ± 1.604 


-28.20 


SDSS J2225-0014 


4.886 


0.005 


Lya forest 


2.144 ±0.885 


-27.39 


SDSS J1616+0501 


4.876 


0.005 


Lya forest 


2.603 ± 1.020 


-27.60 


SDSS J2 147-0838 


4.588 


0.005 


Lya forest 


3.506 ± 1.245 


-27.93 



(Mpc) 



Tbkg 

HT 12 s" 



logfTbkg) 



22.1 
16.8 
22.3 
17.3 
16.3 
13.3 
10 
12 
14.0 
4.0 
8.3 
8.2 
6.0 
3.3 
5.3 



18.2 
10.0 
11.2 
6.7 
+24.0 
11.6 
+ 10.1 

6.4 
+ 13.9 

7.5 
+ 7.9 
5.0 
+5.6 
-3.7 
+6.4 
-4.3 
+8.6 
-5.3 
+10.3 
-2.9 
+ 8.8 
4.3 
+3.0 
2.2 
+4.2 
2.5 
+2.6 
1.5 
+4.0 
2.3 



0.08 



+ 0.19 
0.06 
+ 0.15 



0-08-o.og 

03+° 10 
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0.10 



0.12 
0.23 
0.43 
0.14 
0.20 
1.82 
0.48 
0.95 
1.02 
3 

1.88 



0.06 
+0.30 
0.08 
+ 0.40 
0.15 
+0.65 
0.26 
+ 0.26 
0.09 
+0.40 
0.13 
+23.60 
-1.69 
+ 1.88 
0.38 
+ 1.13 
0.52 
+2.51 
-0.72 
10.80 
2.86 
+4.58 
1.33 



-13.10 ±0.53 
-13.09 ±0.46 
-13.51 ±0.64 
-12.99 ±0.42 
-12.94 ±0.55 
-12.64 ±0.43 
-12.37 ±0.40 
-12.86 ±0.46 
-12.71 ±0.48 
-11.74 ± 1.15 
-12.32 ±0.69 
-12.02 ±0.34 
-11.99 ±0.54 
-11.41 ±0.58 
-11.73 ±0.54 



and high redshift (z > 5.5) samples, and Tbk g was measured by 
finding where Yli AFi — 0, with AF calculated for each spec- 
trum using the same trial Tbkg- For the lower redshift bin, con- 
taining 8 spectra, the average redshift was z = 5.04 and gave 
logfTbkg) = -12.15±0.16 (0.32) at 68 (95) per cent confidence, 
whilst for the higher redshift bin with 7 spectra and an average red- 
shift of z = 6.09, log(r bk g) = -12.84 ± 0.18 (0.36). These will 
subsequently be referred to as the z ~ 5 and z ~ 6 samples, re- 
spectively, for ease of comparison to the simulations and previous 
work. The implications of these results will now be discussed. 



4.2 Comparison to previous work 

The results of fitting multiple lines of sight with a constant UVB 
are plotted in Fig. Q~3] We also plot a selection of UVB estimates 
from the literature, where the literature results have been scaled to 
the cosmology used in this paper, as well as the same temperature - 
density relation parameters used by iBolton & Haehnei 3 J2007bh . 
such that T = T A 7_1 , with T and 7 held constant at 10 4 K 
and 1.3 respectively. The UVB model of lHaardt & MadaJ J200lh 
is also shown. As previously noted in Fig. [TO] whilst the error bars 
are large for the individual sightlines, there is a clear trend of a 
decreasing UVB intensity in the redshift range z ~ 5 — 6. This 
decrease is more pronounced in the binned results. 

The results presented here are th e first proximity e ffect m ea- 
surements of Tbk g at these redshifts. ball' Agli o et al. I d2009h de- 
tected the proximity effect in 1733 spectra from the SDSS over the 
range 2 < z < 4.5 and found that the UVB seemed to be remark- 
ably flat over this redshift range, withlog(r b k g ) = -11.78±0.07. 
This value is consistent with our results from individual sightlines 
over 4.5 < z < 5. A more d i rect c omparison can be made with 
the results o f lDair Aglio et al] (QoOX ). who measured the proxim- 
ity effect in 40 spectra from the UVES instrument on the VLT. 
Comparing our results to the iDair Agli o et al. I d2008l) results may 
be more appropriate, as both studies used high-resolution spectra. 



Combined, the proximity effect results suggest a smooth decrease 
in UVB from z ~ 2 to z ~ 6 by an order of magnitude CFig.l 13b. 

The range z ~ 5 — 6 has previously only been 
probe d using flux decr e ments dMcDonald & Miralda-Escudel 



200 it iMeiksin & White! |2004 IBolton & Haehneltl l2007bT 



Wvithe & Boltonll2010l) . There is excellent agreement between 



binned results and those from flux decrements (see Fig.ll3t, which 
also suggest a significant decline in the UVB from z ~ 4 to z ~ 6. 

Our measured evolution of the UVB intensity has important 
implications for reionization. The H 1 photoionisation rate should 
scale as T(z) oc l(vo, z)e Vo , where 1(vq, z) is the mean free path 
of ionising photons and e„ is the ionising emissivity. The redshift 
evolution of Tbk g therefore gives an insight into the evolution of 
these two key variables. During the 'percolation' stage at the end 
of reionization, l{vo,z) is expected to evolve rapidly, in marked 
contrast to its gradua l evolution in the post-reionization Universe 
dGnedin & Fanl200q) . Consequently, the end of reionization should 
be indicated by a break in the evolution of l(fo, z). The smooth red- 
shift evolution of our measurements of the UVB intensity (see left 
panel of Fig.llOt implies that both l(yo, z) and e„ are also evolv- 
ing smoothly in the redshift range 4.6 < 2 < 6.4, as otherwise 
they would have to both evolve rapidly simultaneously in opposite 
directions, which appears very unlikely. This suggests that percola- 
tion has occu rred at higher redshifts than ar e probed by our sample 
(although see lFurlanetto & Mesinger|[2009l) . 

Our measurements of the UVB can be combined with mea- 
surements of the mean free pat h to place constraints on t he evolu- 
tion of the ionising emissivity. ISongaila & Cowiel d2010h recently 
measured the incidence of LLS over < z < 6. They find an 
evolution in the mean free path which can be approximated as 

. , .,+0.36 

l(u , z) = 50[(1 + z)/4.5]~* . Using this fit, we infer that 

1(vq, z) decreases by a factor ~ 1.5 — 2.5 from z — 5 — 6. At 
face value the decrease of our measurements of Tbkg with redshift 
(a factor of ~ 2.5 — 8) therefore imply an emissivity that is either 
roughly constant or drops by up to a factor of ~ 5. Note that this 
is consistent with the decrease in the (dust-corrected) UV luminos- 
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Figure 12. The 1 5 spectra used in this study, ordered by redshift. The spectra are presented in their rest wavelength, A rcs t , in order to emphasise the proximity 
effect region bluewards of the Lya emission line, on the far right of the plot at 1216 A. They have been normalised with the continuum fitting process described 
in Section |2~T1 and have also been smoothed to a common pixel size of 16.7kms _1 (observed frame) for clarity. The red dashed line indicates the derived 
value of i? cq , whilst the blue dotted line is the expected fall-off in average flux. The area of spectrum to the right of the green dashed line was the section used 
for the proximity effect measurement, and is 40/i~ 1 comoving Mpc long. 
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Figure 13. The evolution of the UVB in the redshift range z = 2 to z = 6. The solid points are the results of this paper, and represent the aggregate analysis 
over multiple sightlines. The spectra were coarsely binned into two subsets, namely z > 5.5 and z < 5.5. The lower redshift subset contained 8 spectra and 
had an average redshift of z av = 5.04 whilst the upper had 7 spectra and 2 av = 6.09. The errors displayed on our points are the 68 per cent confidence 
intervals, and exclude any correction due to the presence of LLS, local overdensities, or a thermal proximity effect. Recent estimates from the literature are 
also plotted with their la errors. The flux decrement results have all been scaled to our adopted cosmology (h, Q m , n b h 2 , erg) = (0.72,0.26,0.024,0.85) 
and to the same tem per ature-density relation, T = TbA 7-1 , with To and 7 held constant at 10 4 K and 1.3 respectively, using the scaling relations from 
iBolton et al] ]2005l) and | Bolton & HaehneltH2007bl) . Our results at both z ~ 5 and 2 ~ 6 agree very well with those from flux decrement measurements. The 
theoretical curve of lHaardt & Madaul feOOln . assuming contributions from quasars and star forming galaxies, is also plotted. 



ity density between 2 = 5 — 6 measured bv lBouwens et al] J2009h 
for rest wavelength ~ 1600 A and integrated to eithe r 0.3L* = a or 
0.04L *^3 (where Lt =3 is the luminosity derived by [Steidel et al] 
dl999h at z ~ 3 and corresponds to Mnoo.AB = -21.07). 

The fit to the mean free path from Songaila & Cowiej d2010h is 
plotted with the measured proximity region sizes in the right panel 
of Fig. [TO] Our proximity r egion sizes become larger than the mean 
free path as measured by ISongaila & Cowiej ( 1201 Of) at z > 5.5. 
If the quasar is able to ionise LLS within the proximity region, 
then the mean free path within the proximity region will increase, 
allowing an enhanced contribution to the ionisation rate from local 
galaxies and AGN. As a result, the proximity region may appear 
larger, and so the UVB may be underestimated in a proximity effect 
analysis. Careful modeling of LLS will be required to determine 
whether this is an important effect. However, it may partly explain 
the decrease in ionising emissivity we infer from z ~ 5 to 6 (see 
also Section [3".4.2| l. 

The ionising emissivity is already very low at 2 = 6 (~ 1.5 
ionising photons per hydrogen atom) and so this evolution can not 
continue to much hi gher redshifts without reion ization failing to 
complete by z = 6 dBolton & H aehnelt 2007b). For this reason, 
unless there is a substantial increase in e„ at very high redshifts, 



the end of reionization, while potentially before z = 6.4, appears 
unlikely to occur much earlier. 



5 CONCLUSIONS 

We have presented new measurements of the photoionisation rate 
of hydrogen by the ultraviolet background (UVB) using the prox- 
imity effect in quasar spectra. The fifteen spectra in the sample 
cover the range 4.6 < z < 6.4, allowing us to conduct the first 
proximity effect measurements of Tbkg at these high redshifts. 

For each quasar, Fbk g was calculated by modelling the total 
photoionisation rate as a function of distance from the quasar, tak- 
ing into account contributions from both the quasar and the UVB. 
The optical depths in the proximity zone were then modified to 
produce the expected optical depths in the absence of the enhanced 
ionisation from the quasar. The preferred value of Tbkg was the one 
for which the resulting mean flux in the proximity zone was equal 
to that of the average Lya forest at the same redshift. 

We investigated a wide range of potential errors and biases 
affecting the proximity effect measurements using full radiative 
transfer simulations. We found the error in Fbkg for an individ- 
ual line of sight to be dominated by the error in the quasar redshift 
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and T c ff . The redshift errors for the sample are generally ^ 0.01, 
with several having a very accurate systemic redshift from CO and 
Mg n . An error in the assumed luminosity of the quasar of up to 40 
per cent would still have a smaller effect than the typical redshift er- 
rors. Our raw measurements of the UVB intensity from individual 
sight lines at z = 5 should be typically overestimated by a factor 
of two, with random errors that are also roughly a factor of two. 
The largest biases result from the finite S/N of the spectra and the 
quasar lying in an overdensity (50 and 30 per cent, respectively). 
At z = 6, the UVB is overestimated by 30 per cent, with a factor 
of 2.5 random error. The largest source of bias at this redshift is the 
uncertainty in T e rf . 

Systematic uncertainties caused by the overdensity of matter 
close to the quasar were found to be smaller at z > 5 than at lower 
redshifts. The effect of the quasars lying in large-scale overdensi- 
ties was found to be small for host halo masses of M ^ 1O 12 2 M , 
the most massive halos in our largest (80/i _1 comoving Mpc) sim- 
ulation. Even in a toy model emulating the density distribution sur- 
rounding a 10 13 Mq halo at z = 6, (and therefore as massive as 
the most massive haloes in the Millennium simulation at this red- 
shift) large-scale overdensities were not dominating the error bud- 
get. This is due to the proximity regions (as defined in this paper) 
being substantially larger than the scale of the overdensities. It is 
however noted that at lower redshifts, where the UVB intensity and 
the mean flux level are higher, proximity regions may well be of 
similar size to the local overdensity. This effect is therefore likely 
very important at z ~ 2 — 3, where discrepancies between mea- 
surem ents from flux decrements and the proximity effect have been 
noted teollinde et al .11 20051 : [Guimaraes et al]|2007t) . 

We also tested the effect of Lyman limit systems (LLS) Using 
the r ecent measurements of th e frequency of LLS at high redshift 
from lSongaila & Cowiej J201C|> . we found that LLS have also only a 
small effect on our analysis, potentially causing us to overestimate 
Tbkg in our combined samples at z ~ 5 and 6 by roughly 0.1 dex. 

Our measurements further assume that the gas temperature 
within the proximity region is similar to that in the general IGM. A 
thermal proximity effect, could have a sizeable impact on our esti- 
mates of the UVB intensity. If the gas in the ~ 5 proper Mpc of the 
proximity region closest to the quasar were e.g. on average 10 4 K 
hotter than the general IGM then the UVB will be underestimated 
by 0.14 (0.21) dex at z — 5 (6). This is comparable in magnitude 
to the biases from other environmental factors, although opposite 
in sign. 

The results quoted in this paper are corrected for these biases, 
except the environmental biases due to a different gas temperature 
and average gas density, and the presence of LLS in the proximity 
zone. We have, however, demonstrated these biases are small or, in 
case of the temperature, are poorly observationally constrained. 

Our measured values of Tbkg, corrected for biases, decline 
significantly from z ~ 5 to 6. For z ~ 5 we find log(rbkg) = 
-12.15 ± 0.16 (0.32) at 68 (95) per cent confidence, whilst at 
z ~ 6 we find log(r bkg ) = -12.77 ± 0.18 (0.36), a decline 
significant at roughly the 3<r level. Within our sample, the UVB 
intensity measured from individual sight lines is seen to decline 
smoothly with redshift over 4.6 < z < 6.4, but show no sign 
of the rapid decline which may be expected in the late stages of 
reionization when there is a rapid change in the attenuation length 
of ionising photons. 

Our results are in good agreement with UVB estimates from 
measurements of the mean flux decrement in the redshift range 
z = 5 — 6, assuming an IGM temperature To = 10 4 K. Both 
the proximity effect and flux decrement measurements imply a de- 



cline in the intensity of the UVB by nearly an order of magnitude 
from z = 4 to z — 6. 

Finally, we have combined our esti mates of Vbks with the evo - 
lution of the mean free path measured bv lSongaila & Cowiej d20ld) . 
At face value the combined measurements imply a decline in the 
ionising emissivity of a factor of about 2 from z ~ 5 to 6, but it is 
important to keep in mind that at these redshifts measurements of 
the mean free path of ionising photons are extremely difficult, and 
that with such a low emissivity reionization could barely have been 
completed by z = 6.4. 

The results presented here represent some of the highest red- 
shift measurements of the UV background made to date, enabling 
us to probe deeper into the late(st) stages of hydrogen reionization. 
The next generation of optical and NIR telescopes will enable ac- 
cess to high signal-to-noise quasar spectra at even higher redshifts, 
leading to improved measurements of the ionising background and 
helping to establish a more complete picture of the final stages of 
the hydrogen reionization epoch. 
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Figure Al. Distribution of the estimated values of log(Tbkg) as a function 
of S/N for 1000 simulated spectra with HIRES resolution. A proximity 
region size of 10 Mpc is assumed. Apart from the luminosity of the quasar 
and the peculiar velocities of the gas, there are no other sources of error. 
The dotted line marks the input value, while the blue dashed line marks the 
mean of the simulated data set. The estimated values of log(Tbkg) are well 
fit by a Gaussian. 



APPENDIX A: DETAILED DISCUSSION OF THE ERROR 
ANALYSIS 

In this Appendix we show how varying some of the parameters 
in Table [3] affects the scatter in log(Fbkg) and induce a system- 
atic bias. By modelling the effects of all these parameters, an ex- 
pected sightline-to-sightline scatter, a exp , and bias, e exp , can be 
calculated for each individual quasar. The former contributes to 
our quoted errors in logfTbkg), and the latter is used such that 



log(rS r g r ) = log(r^) - e exp . All values of log(r bkg ) quoted 
in Table [5] have been corrected for these biases. For most of the 
spectra, the combined correction for all the sources of error is 
fei» 5, 0.2 dex. 



Figure A2. Same as Fig. IA II but for simulated data with MIKE velocity 
resolution. 
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Figure A3. Expected distribution of errors in log(Fbkg) for various errors 
in the luminosity, assuming a fixed S/N of 20, and HIRES resolution. 
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Figure A4. Expected distribution of errors in log(rbkg) for various errors 
in the redshift, for an assumed S/N of 20, and HIRES resolution. The top 
panels assume no redshift error. The next two panels down represent CO 
and Mg II redshifts, respectively, and the bottom panels represent the errors 
on redshifts determined from the onset of the Lyct forest. 



Figure A5. Expected distribution of errors in log(Fbkg) for various input 
proximity region sizes (i.e. luminosity of the quasar), for a fixed S/N of 
20, and HIRES resolution. 
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Figure A6. Expected distribution of errors in log(rbkg) for various quasar 
redshifts, for noiseless spectra (i.e. S/N = 00), and at HIRES resolution. 
Even for noiseless data, the method struggles at lower redshift. The distri- 
bution loses its Gaussian symmetry, although the mode o f the distribution 
does m aintain its position close to the input value (see also lDaH'Aglio etal] 
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